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Abstract 

Electromagnetic emission rates from a thermalized hadronic gas are impor- 
tant for the interpretation of dilepton signals from heavy- ion collisions. Al- 
though there is a consensus in the literature about rates for a pure meson gas, 
qualitative differences appear with a finite baryon density. We show this to be 
essentially due to the way in which the irN background is treated in regards 
to the nucleon resonances. Using a background constrained by unitarity and 
broken chiral symmetry, it is emphasized that the thermalized hadronic gas 
can be considered dilute. 
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Recent relativistic heavy-ion collisions at CERN have reported an excess of dileptons over 
a broad range of invariant mass M = 300-500 MeV Various theoretical assessments of 
the dilepton emission rates have tried to understand this excess. In the absence of baryons, 
the rates achieved using detailed reaction processes ||£| and spectral sum rules |5||| agree 
with each other, but fail to reproduce the enhancement in the data. In the presence of 
baryons, calculations using many-body dynamics J7||| lead to larger rates below the p peak 
compared to results from spectral considerations constrained by broken chiral symmetry || . 

This is illustrated by the left plot of Fig. [1] for pertinent temperature T = 150 MeV and 
baryon chemical potential^ p = 520 MeV. Although neither rate accommodates the data in 
a realistic hydrodynamical evolution [nj, the rates from Ref. [^,|J are about two to three 
times larger than those of Ref. f|, being within only two standard deviations of the data. 
We show below that the discrepancy between these two rates originates chiefly from the ttN 
background, and that the hadronic gas is essentially dilute, in confirmation of our earlier 
work ||. 

Although perturbative unitarity fixes the ttN background uniquely in the Compton am- 
plitude, its extrapolation away from the photon-point requires care. In Ref. |9[] this ex- 
trapolation was implemented by paying due care to threshold unitarity. In the future, the 
complicated issue of background versus resonances should be properly resolved by enlarging 
the theoretical analysis to pion photo-production as well as pion knock-out while obeying 



unitarity and broken chiral symmetry as detailed in Ref. [p |. For now, we concentrate on 
illuminating the reason for the above differences in the two theoretical estimates containing 
nucleons. 

In a thermal equilibrated hadronic gas, the rate R of dileptons produced in an unit four 
volume follows from the thermal expectation value of the electromagnetic current-current 
correlation function [|T^]. For massless leptons with momenta pi and P2, the rate per unit 
invariant momentum q = p± + p2 is given by || 

A = T^m ImWF M (D 

where a = e 2 /4ir is the fine structure constant, and 

W F (q) = % Jd 4 x e iq - x Tr ( e -(H-**N-n)/T T *jm( x ) j^q)) . (2) 

eJ M is the hadronic part of the electromagnetic current, H is the hadronic Hamiltonian, p 
is the baryon chemical potential, N is the baryon number operator, Q is the Gibbs energy, 
T is the temperature, and the trace is over a complete set of hadron states. 

We expand the trace in Eq. (Q) using pion and nucleon states. This is justified for 
temperatures T < m n and final nucleon densities p^ < 3po with p = 0.16 fm -3 , since the 



1 We express our rates using a complete set of stable states § and therefore \i refers solely to 
nucleons. In the analysis of Ref. Pj8|, (i refers to all possible baryons as they retain the unstable 
particles in their final states. The rate comparison done here is only meaningful for fixed T and fi. 
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FIG. 1. The left plot compares nucleon rates from Refs. || (solid) and |Q] (dashed). The right 
plot shows the results for the (isospin averaged) Compton scattering on the nucleon from the A 
(solid) and background (dashed) versus data |jl3| , 

relevant expansion parameter n for each particle is less than | in that regime. This allows 
us to only retain terms up to first order in density |J 

ImW F (g) = -3g 2 Imn y (g 2 ) + -^ f dirW%(q, k) + fdN W%(q,p) + 0(4, k 2 n , k w k n ) . (3) 

The phase space factors are 

d 3 p 1 1 a — d 3 k 1 1 

(2tt) 3 2E p e^v-»)IT + i and d7T ~ (2^p2^ e ^/ T -l' 

with nucleon energy £7 p = \/ m 2 + p 2 and pion energy tUk = ym^ + k 2 . The first term in 
Eq. (0) is the transverse part of the isovector correlator (0|T*VV|0) and summarizes the 
results of the resonance gas model. It is given by e + e~ annihilation data. At low and 
moderate q 2 this is dominated by the p and p' while at high q 2 its tail is determined by the 
qq spectrum. 

The term linear in pion density can be expressed in terms of experimentally measurable 
quantities by use of chiral reduction formulas. The important contributions are [U 

Wf (q, k) ~ 12g 2 Im U v (q 2 ) - 6(k + q) 2 lm U A ((k + qf) + (q -> -q) 

+ 8{(k ■ q) 2 - m\q 2 ) Im Ii v {q 2 ) Re (A R (k + q) + A R (k - q)) (4) 

with Aft(k) the retarded pion propagator and 11^ the transverse part of the iso-axial correla- 
tor (0|T*jyJ,i|0) which follows from tau decay data ||. It is dominated by the ai resonance. 

The term linear in nucleon density is just the spin-averaged forward Compton scattering 
amplitude on the nucleon with virtual photons. This is only measured for various values of 
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q 2 < 0. However, the dilepton and photon rates require q 2 > 0. Therefore, only the photon 
rate for this term can be determined directly from data by use of the optical theorem 



e 



W* N (q,p) = -4( s -m 2 )J2^N(s) (5) 

with s = (p + q) 2 . For off-shell photons, we must resort to chiral constraints to determine 
the nucleon contribution to the dilepton rate. Broken chiral symmetry dictates uniquely the 
form of the strong interaction Lagrangian (at tree level) for spin-| particles. Perturbative 
unitarity follows from an on-shell loop-expansion in 1/ f n that enforces current conservation 
and crossing symmetry. To one-loop, the ttN contribution is parameter free. The large 
contribution of the A to the Compton amplitude near threshold is readily taken into account 
by adding it as a unitarized tree term to the one-loop result M. The result for Compton 
scattering on the nucleon from Ref. |J is shown in the right plot of Fig. [I] versus data ||13|| . 

This fit to the Compton data is good, given the fact it is a parameter-free analysis. 
However, we must determine the role of the iV*(1520), since about 20% of the cross section 
is unaccounted for in that kinematic region. More importantly, it has been suggested that 
the decay of this resonance in matter is what feeds the dilepton rate enhancement [7]|§. 
We note that the iV*(1520) has about a 50% branching ratio to the irN channel, making it 
difficult to disentangle from the irN background. 

The contribution of the iV*(1520) to the Compton amplitude follows readily from the 
transition matrix element 



(N(p)\J IM \N*(k))=u(p) 



Q*(Q 2 ){l^ - Q^i) + iS^q 2 )alq x qu 

+ (R*((i 2 ) + iRM 2 )) (q»Qu - g^q 2 " 



u v {k) . (6) 



The form factors Q*, S 1 *, R* and R* comprise a maximal set and are real by time reversal 
invariance. This decomposition parallels the one adopted for the A in Ref. || with ad- 
justments for the difference in parity and isospin between the two resonances. Hence, the 
iV*(1520) contribution to the Compton amplitude A4*(s, q 2 ) is the same form as Ji4(s, q 2 ) 
for the A quoted in Ref. as long as uia is changed to — m* (parity) and the prefactor | 
is changed to 1 (isospin). Just as for the A, the form factors in Eq. (|6]) will be assumed q 2 
independent. 

The resulting couplings are constrained by the nucleon polarizabilities, the i?/M-ratio, 
and the Compton amplitude. The nucleon polarizabilities can be separated into the loop, 
A, and iV* contributions: 



la m m A + m 2 2a m mi + m 2 

A m 2 A 
la Q 2 2a Q 2 



a + (3 = (a + /3) loop + — — f 2 Q 2 + —— *- 2 Qj , (7) 

9 m\ mi — m z 6 mi mi — m z 



= Aoap + -Q-— _^ — r^L . ( 8 ) 

Note that the A contribution gives E2/M1 ~ —0.3 while empirically it is —0.015, and the 
N* contribution gives E1/M2 ~ —3 while empirically it is between —2 and —3. Using 
that the (isospin averaged) experimental results |13| are ^(15.8 + 14.2) x 10~ 4 fm 3 and 
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FIG. 2. Background contribution (dashed) and total contribution (solid) to the Compton 
amplitude compared to data |13| (crosses). The left plot is for the full background and right plot 
is for one-third the background. 



5(6.0 + 2.1) x 10~ 4 fm 3 and the loop contributions Q are §(8.8 + 5.4) x 10~ 4 fm 3 and 
— §(1.5 + 2.2) x 10~ 4 fm 3 respectively, we may estimate the values of the couplings Q and 
Q*, with some freedom available from error bars. We choose^ Q = 1.8/m and Q* = 2.2/ m 
for which the polarizabilities are a + (3 = 13 x 10~ 4 fm 3 and j3 = 3.6 x 10~ 4 fm 3 . 

We note that i?* and do not contribute at q 2 = 0, and are mainly associated with 
the longitudinal part of the electroproduction cross section for q 2 < 0. By analogy with the 
A, they will be set to zero ||. We are therefore left to constrain S* which is readily done 
by fitting to the Compton amplitude. Given the tcN background, we find = — 1.5/m 2 . 
Similarly S = 1.5/m 2 for the A offers a better fit than the value S = 1.2/ m 2 quoted in 
Ref. |J. The radiative width from the resonant part of the A and N* are 0.27 MeV and 
0.19 MeV respectively, to which the contribution from the irN background should be added. 

The fit to the Compton amplitude resulting from this new set of parameters and including 
the 7rN and A background is shown in the left plot of Fig. 0. The background contribution 
is shown by the dashed line and the total contribution by the solid line. The crosses are 
the data [13|]. The difference with Ref. || is a change in the A couplings (see above) to fit 
the data and the addition of the iV* to account for the 20% discrepancy. The background 
is parameter free. 

To assess the importance of the irN background versus the N* resonance in the electro- 
magnetic emission rates, it is useful to do the following exercise: we drop the background 



2 The value of Q = 2.75 fixed purely from the A as quoted in Ref. || leads to an overestimation 
of the Compton cross section when added to the background, as seen in the right plot of Fig. |l]. 
This comes from a double counting indicative of a treatment without the proper unitarization of 
the background and resonances. 
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FIG. 3. Compton amplitude on nuclei comparing the two background choices to data [jD] ]. 

by a factor of three and change the A and N* parameters to recover the Compton ampli- 
tude, at least through the N* region. This leads to the right plot of Fig. [|, again with 
the background contribution shown by the dashed lines and the total contribution by the 
solid line. The couplings are still fit to the experimental polarizabilities in this case giving: 
(Q,QJ = (1.8,4.3)/m and (5,5,) = (-2.6, -3.0)/m 2 . The E/M-ratios are therefore the 
same as before. This smaller background is consistent with the one used in Ref. ||. This is 
clearly different from our approach in which the background is first constrained by broken 
chiral symmetry. 

The Compton amplitude on nuclei is also measurable |l5j . To leading order in the nucleon 



density it can be readily assessed using the Compton amplitude on a single nucleon smeared 
by Fermi motion (ignoring Pauli blocking), 

with q 2 = and s = m 2 + 2qo(E p — \p\ cos 6^) and pp = 265 MeV. The results from Eq. (^|) 
are shown in Fig. |3] versus the data [[DJ. The solid line follows from the left plot of Fig. |2| 
while the dashed line from the right plot of Fig. |2|. Since we do not account for particle-hole 
excitations below the ttN threshold, we have not reproduced the large increase near zero 
go- This has no consequence for the dilepton rates to be discussed below. The discrepancy 
of the dashed curve at high go follows from neglecting the higher resonances in Fig. 
Fermi motion smears the contribution of the A and iV*(1520) resonances more than the 
background. The results for the smaller background (dashed line) is in better agreement 
with the data, although both predictions are within one standard deviation of the data. 
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FIG. 4. The left shows the dilepton rate for T = 150 MeV and p = 540 MeV comparing our 
old result with that of the one-third background and full background. The right shows the ratio 
of the new results to the old result. 

The nucleon contribution to the Compton amplitude as shown in Fig. can be used in 
Eq. (|3|) in conjunction with the pion contribution as discussed in Ref. M. Using the identity 
(M = V?) 

dR _ 2 dR 

~¥q ~ ndM 2 dydq 2 ± ' 

and integrating over rapidity y and transverse momentum qj_, the dilepton rates are shown 
in the left plot of Fig. |] for T = 150 MeV and p = 540 MeV (corresponding to a nucleon 
density of po). The purely pionic contribution is shown by the dotted line, the dot-dashed 
line is the result corresponding to the left plot of Fig. while the dashed line is the result 
corresponding to the right plot of Fig. ^. The result quoted in Ref. |J is shown by the 
solid line. As expected, the changes between the solid and dot-dashed lines are small and 
within expectations. More importantly, the dashed rate, following from the smaller back- 
ground contribution, is substantially larger due to the enhancement caused by the additional 
strength of the N* resonance^. 

To make further comparison between the dilepton rates from different backgrounds, the 
right plot of Fig. ^ shows the ratios of the new rates to the old rate of Ref. ||. The dot- 
dashed line is for the rate associated with the left plot of Fig. |2] while the dashed line is 
for the rate associated with the right plot of Fig. |]. It shows that in the intermediate 
mass range, the enhancement of the rates obtained from a smaller background || can be 



3 There is also an additional strength from the A resonance, but it is significantly less. 



7 



two times more than those from a larger^ background 0. The photon emission rates show 
similar behavior to the dilepton rates, with the smaller background giving twice as many 
photons as the larger one. 

Medium modifications as discussed in Ref. J7]|| (and by others [fL6l ) in addition appear 
to deplete the dilepton rates in the p region, which is qualitatively different from our rates 
with any background. The iV* contribution is fully within the theoretical uncertainties for 
our old rates as long as the full background is used. This implies that in the low mass 
range, the rates are well accounted for by merely a dilute gas. The discrepancy between our 
rates and those in Ref. [[7],[5| around the p is real and may be indicative of some coherence 
or lack thereof. This point can be sorted out through higher resolution measurements as we 
suggested in Ref. |5]||. 

The electromagnetic emission rates in a baryon-free gas are well established theoreti- 
cally HH . In a baryon-rich gas, the discrepancy between the rates reported in Ref. and 
those in Ref. || is mainly due to the choice of the irN background. The Compton amplitude 
on the nucleon and nuclei can both accommodate either a small || or large || ttN back- 
ground. A calculation with a small background requires more strength in the resonances, 
leading to dilepton and photon rates that are about two times larger than from a large 
background at T ~ m n and p^ ~ Po- 

Perturbative unitarity and broken chiral symmetry uniquely specify the larger ttN back- 
ground. A one-loop expression for this quantity was used here. Adding more loops can 
alter and possibly deplete the background above the two-pion threshold. However, these 
corrections are chirally suppressed as evidenced by the fact that the one-loop calculation is 
a good approximation to the background above the resonances. A systematic understanding 
of how the higher order contributions enter, like through pion photon-production and pion 
knock-out, is still important. This can be done reliably using the framework proposed in 



Ref. ||11|| , where due care is paid to unitarity and broken chiral symmetry. As always, it 
is important to constrain the dynamics of relativistic heavy ion collisions by broken chiral 
symmetry so that novel effects can be reliably assessed. 
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4 We have checked that an adjustment of the background to coincide with Ref. |7[], which is about 
75% of that used in Ref. || , yields a ratio comparable to the full background as shown in the right 
plot of Fig. @ 
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